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FOTO PORTADA 





San Carlos de Bariloche - Argentina 


Imagen de satélite QuickBird de 
DigitalGlobe, producto Pansharpened 
resolución espacial 0.60 metros, en 
composición Color Natural (3,2,1) de la 
Ciudad de San Carlos de Bariloche. Esta 
majestuosa ciudad se encuentra a orillas del 
lago Nahuel Huapi, espejo de agua que junto 
a montaña y bosques integra el contexto de 
uno de los lugares más bellos del mundo. 
Esta ubicada al oeste de la provincia de Río 
Negro y a 1.640 kilómetros de la ciudad de 
Buenos Aires. 

En la imagen se observa el Centro Cívico de 
Bariloche y el puerto de la ciudad a orillas del 
Lago Nahuel Huapi. 

Esta ciudad fue sede del VW Simposio 
Latinoamericano de Percepción Remota, en 
elaño 1989. 
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PLAN EDITORIAL SELPER 


A partir de las decisiones adoptadas en el marco del X 
Simposio Latinoamericano de Percepción Remota y Sistemas de 
Información Espacial, organizado por nuestra Sociedad; llevado 
a cabo en Cochabamba, Bolivia, en noviembre de 2002, el 
Comité Editorial ha tomado el compromiso de regularizar las 
publicaciones de la Revista Técnica Científica. 


En este sentido es necesario destacar el invalorable apoyo 
obtenido de la Agencia Espacial Europea (ESA) por la impresión 
de los volúmenes correspondientes a los años 2000, 2001, 2002, 
2003, 2004, 2005 y 2006 respectivamente, los cuales se 
encuentran disponibles para los miembros de nuestra Sociedad. 


La 1% edición del año 2007 se llevó a cabo gracias al importante 
aporte recibido por parte del Capítulo Nacional Selper México. 


En esta oportunidad hacemos llegar la publicación del volumen 
26 número 1, donde se incluyen trabajos que han sido enviados 
por especialistas en el área de Teledetección y SIG que 
desarrollan sus actividades en diferentes instituciones de 
América Latina, los cuales han sido remitidos a expertos 
evaluadores para su revisión y aprobación. 


Nos toca agradecer muy especialmente la labor llevada a cabo 
por el grupo de evaluadores que ha participado en la revisión de 


los mismos. 


NORMAS PARA LOS AUTORES 


Los artículos recibidos serán enviados a tres (3) 
expertos en la temática para su revisión. Los 
trabajos aprobados serán publicados en estricto 
orden, de acuerdo a las fechas de llegada de las 
contribuciones. 

Los idiomas oficiales SELPER son: Español, 
Portugués e Inglés. 

Los trabajos deberán estructurarse 
contemplando las siguientes secciones: 


a) Título del trabajo. Nombre de los autores y 
direcciones completas 

b) Resumen (no más de 150 palabras) 
indicando al final las palabras claves. Deberá 
incluirse en Español o Portugués, además de 
Inglés 

c) Introducción 

d) Objetivos 

e) Metodología empleada y materiales 


f) Resultados obtenidos 

g) Conclusiones 

h) Bibliografía: sólo se incluirá la citada en el 
texto. Se indicarán los autores, por orden 
alfabético, año, título, revista o medio donde fue 
publicado, incluyendo volumen y páginas, cuando 
corresponda. 


Los títulos y subtítulos de cada sección deberán 
estar claramente indicados (ya sea con 
numeración o tamaño de letras). Las tablas, fotos 
y figuras deberán ser suficientemente nítidas, 
llevar un título y estar numeradas en forma 
consecutiva. 

Se deberá enviar una copia del trabajo en 
formato Word y una copia papel. La extensión 
total del trabajo no deberá superar las 12 páginas, 
(DIN-A4). 

Los trabajos se enviarán a: editorial(Wselper.org 
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RESUMEN 


La evapotranspiración (ET) constituye el segundo 
elemento del ciclo hidrológico del agua en la 
atmósfera, a ella están ligados diversos y complejos 
procesos en el que intervienen agua, suelo planta así 
como otros aspectos medioambientales y de manejo, 
su determinación y utilización hace imperante por 
especialistas que planifican, toman decisiones, 
explotan los recursos hídricos, ambientales y las 
geociencias en general. 

En los años 90 la FAO establece el método de 
Penman-Monteith como procedimiento estándar para 
su determinación donde se requiere de observa- 
ciones sinópticas de estaciones meteoro-lógicas las 
cuales generalmente no están distribuidas 
espacialmente que mediante herramientas 
geoestadísticas se pudiera suplir de información 
precisa de algunas zonas adyacentes con un nivel de 
confianza requerido. 

En Cuba uno de lo trabajos mas relevantes y actuales 
son los desarrollados por el Instituto de Meteorología 
(Solano 2003), donde publica las tasas de 
evapotranspiración de referencia del Territorio 
Nacional, el presente trabajo pretende complementar 
estos estudios ofreciendo una metodología de trabajo 
para la obtención de mapas de intervalos de 
confianza dotando al estudio de resultados 
cualitativos. 


INTRODUCCIÓN 


La evapotranspiración (ET) constituye el segundo 
elemento del ciclo hidrológico del agua en la 
atmósfera, a ella están ligados diversos y complejos 

procesos en el que intervienen agua, suelo planta así 
como otros aspectos medioambientales y de manejo; 
su determinación y uso se hace imperante 
porespecialistas que planifican, toman decisiones y 
explotan los recursos hídricos, ambientales y las 
geociencias en general. Desde siglos pasados se han 
generado una serie de modelos biofísicos y empíricos 
para lograr su aproximación pero siempre ha existido 
una limitación que es la alta variabilidad espacial y 
temporal de cada uno de sus componentes, en los 


años 90 la FAO establece el método de Penman- 
Monteith como procedimientoestándar para su 
determinación donde se requiere de observaciones 
sinópticas de estaciones meteorológicas, las cuales 
no siempre están disponibles ni distribuidas 
espacialmente de forma que mediante 
interpolaciones se pudiera suplir de información de 
algunas zonas adyacentes, con un nivel de confianza 
requerido o al menos informarlo. 

El área de Cuba es de110 860 Km2 (ACC 1989), de 
acuerdo a informaciones del Instituto de Meteorología 
existen unas 68 estaciones climáticas de diversas 
categorías que realizan observaciones del clima de 
ellas 62 pueden proveer de informaciones necesarias 
para la estimación de de la evapotranspiración de 
referencia lo que daría una densidad de cubrimiento 
de 1788 K2, estas estaciones fueron situadas no 
atendiendo a un objetivo agroclimático (algunas de 
ellas situadas en puntos como cayos, costas etc.) por 
tal motivo si se pretende realizar un análisis para 
áreas de interés agropecuario debe tenerse en 
cuenta cuales deben seleccionarse a fin de no 
generalizar información (propagación de errores), al 
realizar una interpolación utilizando cualquiera de los 
métodos existentes siempre el número de datos 
disponibles y su distribución espacial constituye la 
primera limitante, pero haciendo un análisis objetivo 
para poder realizar cualquier estudio relacionado con 
la evapotranspiración de referencia es obligado a 
recurrir a esta fuente de información al no existir otra, 
si se consulta información a cerca de la variabilidad 
espacial y temporal de la evapotranspiración en Cuba, 
las publicaciones son insuficientes y algunas carecen 
de un nivel de robustez requerido ya que están 
basados en análisis de datos los cuales fueron 
generalizados a niveles mensuales y de períodos 
hidrológicos donde se enmascara la verdadera 
naturaleza de la variabilidad espacial y temporal de la 
ET; la FAO recomienda la utilización de informaciones 
decenales para este tipo de estudio (Smith 1990), 
entre los trabajos relevantes son de obligada consulta 
los desarrollados por el Instituto de Meteorología 
"Zonificación de la evapotranspiración de referencia 
en Cuba. Atlas Agrometeorológico de 
disponibilidades hídricas para una agricultura de 
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secano" (Solano 2003) al no existir otros 
antecedentes basados en análisis de series a nivel 
decenal en el país. 

El presente trabajo tiene como objetivo 
complementar este estudio estableciendo una 
metodología para el manejo de datos geoespaciales 
de esta naturaleza haciendo un análisis espacio 
temporal de esta variable y determinando los 
intervalos de confianza para las tasas de 
evapotranspiración de referencia, dotando al mismo 
de un análisis cualitativo necesario para la toma de 
decisiones económicas, ambientales y de seguridad 
alimentaría. 


MATERIALES Y MÉTODOS 


Para el presente estudio fueron utilizados datos 
decenales de series históricas (30 años), mapas e 
informaciones obtenidas del Atlas Agrometeorológico 
de disponibilidades hídricas para una agricultura de 
secano publicados por el Instituto de Meteorología 
(Solano 2003) (cortesía de sus autores), donde se 
ofrecen los datos de ETo para el Territorio Nacional de 
62 estaciones climáticas determinados por el 
procedimiento estándar de FAO-Penman-Monteith 
(Allen 1998). 





FIG. 1. Localización de las estaciones 
climáticas en Cuba 


Para la aplicación de la metodología es necesario 
realizar una estratificación de los valores de ETo 
mediante la geomorfología del relieve por la influencia 
que este aporta a la variable de estudio (Smith 1990), 
por este motivo se necesitan informaciones de esta 
naturaleza. Al consultar información al respecto de 
los estudios realizados en Cuba existen diferentes 
criterios (Hernández S.J 1993) señala que los 
resultados actuales de investigación de la 
geomorfología cubana obtenidos hasta hoy, reflejan 
un grado desigual e insufic iente de estudio dado por 
diversas causas. Uno de los trabajos más recientes 
que considera los estudios antes realizados y utiliza 
técnicas más avanzadas es el desarrollado por Reyes 
(2004) donde utiliza procesamientos utilizando 
técnicas de radar (Shuttle Radar Topography Misión. 
Jet Propulsión ) con resoluciones espaciales de 90 m 
(NASA 2002). 

El modelo digital de elevaciones propuesto por su 
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autor fue el seleccionado y utilizado (cortesía de su 
autor) en los posteriores análisis de este trabajo. 
Como procedimiento de trabajo se estableció el 
siguiente: 


: Análisis estadístico exploratorio descriptivo de la 
totalidad de los datos sin realizar ninguna 
transformación. 


: Prueba de Normalidad de Anderson - Darling la cual 
es recomendada como la prueba más robusta 
existente y valida cuando se quiere precisar los 
limites de la serie en cuanto a la normalidad, varios 
autores (Brown 1974; Filliben 1975; Sheather 1986) 
la consideran superior a las pruebas de Ryan-Joiner y 
Kolmogorov-Smirnov. 


: Análisis grafo - analítico de la distribución espacial y 
temporal de la ETo. 


- Estratificación de las poblaciones de la series de ETo 
a partirde la geomorfología . 


: Establecimiento de las nuevas clases a partir de las 
series de ETo obtenidas. 


- Verificación de la normalidad de las clases. 


«Análisis de patrones para identificar el 
comportamiento espacial de los datos. 


- Cálculo de los semi-variogramas experimentales. 


Modelación de los semi-variogramas y 
establecimiento de los modelos. 


: Interpolación mediante el método del Kriging. 
- Cálculos de los errores (desviación estándar). 


: Estadística descriptiva y prueba de normalidad de 
los errores. 


- Establecimiento de los intervalos de confianza. 


-Mapificación y despliegue de los resultados 
mediante sistemas de información geográficos (SIG). 


ANÁLISIS Y DISCUSIÓN DE LOS RESULTADOS 


Estadística exploratoria descriptiva: Se realizó un 
análisis estadístico exploratorio con la totalidad de los 
datos sin realizar ninguna transformación según lo 
que sugiere Rossiter (2003) para determinar la 
naturaleza de los datos y en qué sentido se debe 
orientar los siguientes pasos de la investigación . 
(Rossiter 2003) La Prueba de Normalidad de 
Anderson - Darling en todos los casos muestra 
suficiente evidencia para afirmar que ninguna de las 
series se comportan como una distribución normal 
para el nivel de significación (a=0.05), entre las 
muestras se observa que existen diferencias entre las 
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Fig.2 Histograma ETo de todas las series 
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Fig. 3 Evapotranspiración de refe rencia por estaciones 
climáticas 


Tabla 1. Estadística descriptiva exploratoria de todas las series 
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medias (Min.= 26.871 y Máx.=52.853) con 
coeficientes de variación (CV min.= 6.590 y CV Máx.= 
11.810) provocadas a que los datos se distribuyen en 
diferentes poblaciones como se muestra en el 
histograma de la Fig. 2 lo que permite suponer que los 
datos obedecen a una distribución que pudiera 
depender de la va riabilidad espacial (X, Y, Z) y 
temporal (T) como se muestra en la Fig. 3, este 
análisis define 3 direcciones de la investigación: 

1-Los análisis deben considerar la variabilidad 
temporal de la evapotranspiración de referencia para 
lo cual hay que decidir cual es el periodo en el cual se 
establecen los valores más representativos de la 
serie . 

2-Los datos pueden estar influenciados por una 
variabilidad espacial (X, Y, Z) que pudiera determinar 
su agrupamiento. 

3- Las series no presentan un comportamiento 
Normal lo que implica su transformación para poder 
emplear métodos de interpolación basados en el 
supuesto de obedecer a una distribución Gausiana. 
Para decidir el período donde la variabilidad temporal 
sea representativa de la serie, al analizar la Fig.3 se 
muestra la alta variabilidad temporal de la ET aún 
cuando se utilizaran series diarias u horarias. 

Sobre la variabilidad temporal de la evapo- 
transpiración existen pocos trabajos realizados, 
sobre esto en un panel de expertos agrupados por la 
FAO en Roma, Smith (1990) señalaba que los 
procedimientos para analizar la varianza estadística 
de las precipitaciones se ha convertido en calculo de 
rutina, tales procedimientos no han sido llevados a 
cabo en el caso de la evapotranspiración y otras 
variables climáticas de un comportamiento no normal, 
estudios comparativos adicionales para estimar la 
variabilidad temporal y recomendar los 
procedimientos prácticos en la planificación y 
explotación de los recursos hídricos necesitan ser 
realizados. 

Este panel de expertos establecía tiempos escala en 
el estudio de la evapotranspiración y señalaba que 
estos están en dependencia del propósito de la 
aplicación que se quiera realizar y la precisión del 
método de predicción con que se trabaje 
recomendando intervalos diarios, decenales y 
mensuales para el método utilizado.(Smith 1990) 

Por lo antes señalado y teniendo en cuenta que el 
presente estudio pretende obtención de valores de 
ET con fines de planeamiento hidráulico y explotación 
para usos de regadíos estableceremos el tiempo 
escala de 10 días. 


Variabilidad espacial de la ET 


El panel de expertos antes mencionado señalaba a su 
vez que el uso de datos climáticos de estaciones 
meteorológicas para áreas agrícolas presenta 
diferentes limitaciones, los procedimientos en que se 
trabaje para la determinación de la evapo- 
transpiración deben tomar en cuenta las diferencias 
micro climáticas, vegetación, suelos y cambios de 
altitud. 
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La introducción del Lapso Adiabático para varias 
alturas fue introducido por Hargreaves (1980) y el 
factor de índice de aridez presentado por Allen y Pruitt 
(1986) son útiles ejemplos de esto; el panel señalaba 
además que las técnicas de teledetección pueden 
proveer de adicional información sobre la variabilidad 
espacial (Smith 1990). Por este motivo la 
estratificación de las series de ETo se realizó teniendo 
en cuenta la geomorfología del relieve atendiendo a la 
clasificación antes mencionada (Reyes 2004). 

Sobre la transformación de series para su 
normalización algunos autores (Rossiter 2003) 
sugieren trabajar con datos sin transformar 
estratificándolos de forma tal que los agrupe en 
poblaciones cuyos extremos se encuentren dentro de 
una serie normalizada, respondiendo a un principio 
biofísico y no matemático. 

Los datos de evapotranspiración de referencia fueron 
estratificados quedando normalizados excepto en el 
primer estrato (0 — 40 m) donde fue necesario 
reclasificarlo en 2 partes (0 -15 m) y (16 - 40 m) ya que 
aun existen 2 poblaciones en una misma serie, en el 
caso de las series 85-125, 459-800 y 800-1500 se 
agruparon normalmente dentro de una misma serie. 
Las diferentes series quedaron agrupadas en las 
clases que se muestran en la Tabla No 2. 


Tabla 2. Clases obtenidas de la estratificación 
geomorfológica 






Mo. de 
datos | Clases 


: No. de 
dos 
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156-1500 


Para la estimación de la evapotranspiración (ETo) la 
determinación del número de datos necesarios para 
que los resultados se encuentren dentro del intervalo 
de confianza (en este caso un 95 % de probabilidad), 
lo deseable es que la ETo se estime con un máximo 
valor de confiabilidad y teniendo en cuenta que esta 
variable está expresada en lámina (1 mm-— 10 m3/ha) 
consideraremos 0.5 mm de margen de error. 


Intervalo de confianza para la población de una media 


niz* << 
td (1) 


Donde: 
Valores medios de la serie 
Valores tabulados de la distribución normal para un 
nivel de confianza (1.960 para un 95 % de 
probabilidad) 
7: Desviación estándar de la serie 
n: Número de observaciones 
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Como se observa en la tabla No. 3 en todas las clases 
es necesario un número de datos superior a los 
existentes para obtener un margen de error adecuado 
lo que implicaría que existieran en el país unas 684 
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estaciones climáticas con informaciones para estimar 
la evapotranspiración de referencia, estas deberían 
estar situadas a su vez con una distribución espacial 
que se lograra una densidad adecuada en cada lugar. 


Tabla 3 Número de datos necesarios para un margen de error de 0.5 mm de 
ETo de intervalo de confianza de la media para un 95 % de probabilidad 


Desv. intervalos de No. de datos 
Clases estándar confianza ' Disponibles | Necesarios 
015 41.654 2 969 43 533 39 774 14 
OCT IRON ECO IET II ECT 
41-85 39671 | 1927 | aoss7 | 38496 | 18 | 68 — 
86-1500 37.708 
Total 


Trabajar con datos insuficientes puede provocar 
márgenes de errores (Tabla No.4) con valores 
sobreestimados o subestimados de evapotrans- 
piración. Por otra parte si se quiere lograr una 
interpolación precisa utilizando técnicas 
geoestadisticas como el Kriging (interpolador óptimo) 
sería necesario estimar de forma fiable los 
semivariogramas experimentales lo que requiere un 
tamaño de muestras mínimo. Gascuel (1994) en su 
trabajo sobre la variabilidad de estimaciones 
espaciales considera 150 datos y apunta como un 
número suficiente, mientras Webster (1992) 
considera esta cantidad como el límite inferior y 
recomienda un tamaño de muestra de 225. 


Tabla 4. Cuantificación en unidades de norma de riego del 
margen de error 


Margen de error 











+ 0.99 + 9.93 
| s099 | «998 | 


Análisis de patrones 


En las interpolaciones se asume que existe una 
aleatoriedad espacial de los puntos, para examinar el 
comportamiento de los datos de puntos en el espacio 
e investigar si estos están distribuidos al azar se 
utiliza el análisis de patrones (Fig 4), dos técnicas 
básicas se utilizaron para examinar el modelo 
espacial de puntos, distancia de vecinos más 
cercanos y reflexivos vecinos más cercanos. 








Fig. 4. Análisis de patrones de la evapotranspiración de 
referencia. 


En las figuras (Fig.5-9) cual se muestran los 
semivariogramas experimentales para cada clase 
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Fig. 5 0-15 m 
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De acuerdo a los semivariogramas se observan que 
existen dependencias espaciales entre los rangos de 
0-400 Km. respecto a la evapotranspiración de 
referencia, esta información puede servir de base en 
la distribución de las estaciones climáticas, por otra 
parte la semivarianza oscila entre rangos de 4-12.5 
mm que implican desviaciones estándar de la media 


ANÁLISIS ESPACIO-TEMPORAL DE LA EVAPOTRANSPIRACIÓN... 





Fig.8 85-1500 m 


Pr ——Z-— 


REVISTA 


SELPER 





e - o » osa e a e ass e e ps te 


Fig. 7 41-85 m 





"-————Vocelo Esferico Enera 1 Clase 9.15 
Ni tel Esglerco Ensro + Cass 15-40 
————Nocelo Esferco Ensra 1 Clase 41-95 p mi 
———_Mecelo Esterica Enero + Clase 351500 





| 


ed 








mn Vai 





























Fig. 9 Semivariograma experimental todas las clases. 


2- 3.5 mm altos si se compara con los valores medios 
que se sitúan entre Min.= 2.6871 mm y Máx.=5.2853 
mm en este rango de distancias. 

En las Fig. 10-19 se muestran los resultados 
obtenidos de la interpolación de la evapotrans- 
piración de referencia mediante un Kriging ordinario y 
sus respectivos mapas de errores. 


Tabla 5 Parámetros de los semivariograma para todas las clases 
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Fig. 16-17 Clase 85-1500m 
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Fig. 18-19 Evapotranspiración de referencia I decena mes de enero y mapa de errores. 


A partir de la combinación de los mapas estimados comportamiento normal de los errores, los valores 
mediante el Kriging ordinario y sus mapas de errores críticos de intervalos de confianza se obtuvieron de la 
se crearon mapas de intervalos de confianza tabla de valores estándar de la curva norma (Tabla 6)!. 


mediante mapas de cálculos, verificando el 
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Tabla 6 Intervalos de confianza y valores críticos de la curva normal. 





4 ICc*rOo 

Donde ¿£¿ son los valores de Kriging estimado y (7 el 
error estimado (estándar desviación). La 
multiplicación del factor c (Valor critico) por los valores 
estimados en el mapa de errores ( de una cara de los 
intervalos de confianza provoca la obtención de los 
mapas de intervalos de confianza para la media. 

Del histograma del mapa de errores se evalúa la 
estadística descriptiva y la prueba de Anderson — 
Darling para verificar la normalidad, en los gráficos de 
las Fig., 20 y Fig.21 Se muestra la normalidad de los 
errores lo que permite estimar los intervalos de Fig. 22 Intervalos de confianza de la evapotranspiración 
confianza. de referencia 





Er tall a a Dr rt rs maga de errors A . 
Los resultados del análisis del mapa de intervalos de 


confianza y el histograma (Fig 22 y 23), nos muestran 
zonas de "incertidumbres " donde los intervalos de 
confianza presentan desviaciones entre 5-6 mm, esto 
coincide con las zonas donde la densidad de 
estaciones climáticas es baja e inclusive nula 
situadas en las zonas oeste de la Isla en las regiones 
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Fig. 20 Estadística descriptiva mapa de Errores regiones del mismo comportamiento es la región mas 


oriental de la isla en la Provincia de Guantánamo en 
los municipios de Yateras, San Antonio del Sur, 
Manuel Tames, Baracoa, Imias y Maisi todas estas 
áreas representan un 7.2 % del territorio nacional, en 
las restantes zonas del país unas 10,713,750 ha se 
comportan dentro de los rangos de intervalos de 
confianza de 0.5-4 mm, y 

representan el 92.6 % de las áreas de Cuba. 
Diferentes autores coinciden que las medidas de 
teledetección espacial permiten extender los 
modelos de cálculo a zonas amplias donde no se 
dispone de datos meteorológicos (Rivas 2004; Rivas. 
R 2004) y proponen una serie de modelos basados en 
el estudio de diferentes sensores y variables 
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Fig. 21 Intervalos de confianza de la Prueba Anderson 
Darling 
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En la figura 22 se muestra el mapa de los intervalos de 
confianza para la media de la evapotranspiración de 
referencia de la primera decena de enero evaluado 
para un 95% de probabilidad. 


Fig. 23 Histograma Intervalos de confianza 


1 3 VOL. 26 N* 1 
JUNIO 2008 


obtenidas a partir de técnicas de teledetección (Jiang 
2003; Batra 2005; Batra 

2006; Wang 2006; Bisht 2005). 

Por otra parte existe una tendencia mundial a la 
utilización sistemas inteligentes para la modelación 
de la evapotranspiración debido a la imposibilidad de 
medir directamente y de forma precisa (error 
instrumental), este parámetro, así como limitaciones 
en cuanto al conocimiento de las leyes biofísicas que 
rigen este comportamiento, sobre este aspecto 
Obregón (2002) (2003) señalaba que el estudio de 
sistemas complejos como es el caso de la 
evapotranspiración tal como se puede evidenciar no 
sólo en la “textura” de las diferentes variables de 
entrada y salida, sino también en el carácter no lineal 
de las funciones de transferencia y heterogeneidad 
de sus parámetros, los modelos empleados en 
ingeniería procuran incorporar expresiones que 
consideran la complejidad de dichos sistemas. Para 
tal efecto se acostumbra a generar aproximaciones 
que pueden ser empleadas en los diferentes 
ejercicios de modelación. Una primera discriminación 
incluye los denominados modelos de “caja negra” y 
los de “caja blanca”. En los primeros la característica 
principal es la ausencia de algún principio físico como 
por ejemplo alguna de las leyes de conservación 
(masa, momento o energía). Para los segundos por lo 
menos se incluye una o varias de estas leyes. En los 
primeros es común también hablar de dos tipos de 
modelos aquellos que se encuentran guiados por 
datos y los basados en el conocimiento. Como 
herramientas de éstos se destacan los denominados 
SISTEMAS INTELIGENTES. 


CONCLUSIONES Y RECOMENDACIONES 


A la vista de los resultados obtenidos podemos 
concluir que:En las condiciones de Cuba la 
evapotranspiración de referencia presenta una alta 
variabilidad espacial (x,y,z) y temporal (t) coincidente 
con lo reportado por la bibliografía internacional, lo 


REFERENCIAS 


ACC, Ed. (1989). Nuevo Atlas Nacional de Cuba, 
Editado por el Instituto de Geografía de la Academia 
de Ciencias de Cuba, el Instituto Cubano de 
Geodesia y Cartografía y el Instituto Geográfico 
Nacional de España 


Allen, R., Pereira, L., Raes, D. and Smith, M. (1998). 
Evapotranspiración del cultivo. Guía para la 
determinación de los requerimientos de agua de los 
cultivos. Estudio FAO Riego y Drenaje No. 56. FAO: 
298. 


Batra, N. (2005). Estimation and comparation of 
Evapotranspiration from multiple satellites for clear 
sky days. Department of Civil 8 Environmental 
Engineering ofthe College of Engineering. Ohio, USA. , 


ANÁLISIS ESPACIO-TEMPORAL DE LA EVAPOTRANSPIRACIÓN... 


REVISTA 


SELPER 


que limita el uso de datos puntuales de estaciones 
climáticas cuando estas no se encuentran con una 
densidad adecuada y distribuidas uniformemente, 
aún cuando se utilicen procedimientos y métodos 
geoestad?sticos avanzados (interpolador óptimo) lo 
que constituye una limitante al sobrestimar o 
subestimar los valores reales de la misma, debe 
restarse atención su utilización para evitar la 
propagación de errores o al menos informar el grado 
de confiabilidad con que se trabaja. 

1. Es importante considerar las escalas de tiempo en 
la estimación de la evapotranspiración como lo 
recomienda la literatura internacional y no trabajar en 
escalas superiores de una década, la generalización 
de información a otras escalas utilizando estadígrafos 
de tendencia central como la media, desviación 
estándar y pruebas de normalidad enmascaran el 
verdadero comportamiento de la serie. 

2. La metodología ofrecida puede servir de herra- 
mienta en el análisis y procesamiento de este tipo de 
información que presenta un carácter estocástico, 
probabilístico y difuso, la utilización de otros 
procedimientos estándar pueden inducirle otras 
fuentes de error. 

3. Este estudio no es conclusivo solo tiene un carácter 
metodológico encaminado a proponer este 
procedimiento de trabajo para la evaluación posterior 
de todas las series históricas de la base de datos con 
que se cuenta pero es indicativo de similares 
resultados a obtener al analizar otros períodos ya que 
los intervalos de confianza están determinados por la 
cantidad de estaciones climáticas existentes, su 
distribución geográfica y su densidad. 

4. Recomendamos la combinación de estos 
procedimientos de trabajo con técnicas de 
teledetección espacial modelados a partir de 
sistemas inteligentes para tratar de reducir las 
incertidumbres de algunas zonas donde la estimación 
no ha alcanzado los valores deseados. 
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ABSTRACT 


Object-based image analysis classifies image objects 
which are created by image segmentation. Image 
segmentation delineates image objects according to a 
certain homogeneity criteria, and at the same time, 
requires spatial contingency. Since objects are 
comprised of a group of pixels, the objects properties 
such as mean, standard deviation, shape, texture, etc. 
can thus be calculated and used in the image analysis. 
Object-based image analysis uses object's spectral 
and spatial information, and information pertaining to 
an object's sub-or super-objects if a multi-level image 
object hierarchy has been created. With all these 
available features, itis critical the proper 

selection of features and thresholds to differentiate 
the classes of interest. Normally, this selection of the 
proper features and thresholds relies on the expertise 
of the analyst and is performed on a trial and error 
process. In this paper, a SEaTH algorithm was used to 
calculate the proper features and thresholds based on 
the properties of the selected representative objects 
of class. This paper performed land cover 
classification using object based image analysis on 
Landsat-7 ETM+ imagery in a mountainous area with 
complex land cover classes. In the definition of the 
classifiers ofthe land cover classes, the used features 
and thresholds were calculated by SEaTH algorithm. 
Eight land cover classes were produced: temperate 
forest, tropical dry forest, orchards, grassland, 
irrigated agriculture, rain-fed agriculture, lava flow, 
and human settleements. In this paper, it is 
difterentiated between object primitives and object of 
interest. Object primitives are building blocks in object 
based image analysis and are produced by image 
segmentation; object of interest here corresponds to 
the earth objects. 

Key words: Multi-resolution segmentation, Object- 
based image analysis, Separability and thresholds. 


1. INTRODUCTION 
Pixel-based image classification usually produces 


highly freckled results in classifying images of 
heterogeneous areas. The resulting thematic maps 
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bear either very general land-cover information, or 
else detailed information with limited accuracies 
(Franklin et al. 2003). In object based image analysis, 
homogeneous image objects at a chosen resolution 
are first extracted and subsequently classified and the 
results produced are smoother and closer to the 
human vision. The development of object based 
image analysis stems primarily from the desire to use 
the semantic information to interpret an image, which 
is not presented in single pixels but rather in 
meaningful objects and their mutual relations. Objects 
are composed of groups of homogeneous pixels, and 
thus more information than just the grey value can be 
used in the classification, which considerably extends 
the possibilities for image analysis. Object based 
image analysis intends to do what an image 
interpreter does: recognizing the colour, shapes, 
textures and coherent regions presented within an 
image and associating meaningful objects and their 
contextual relations. 

A feature represents certain information concerning 
objects of interest, such as spectral related 
information, shape, texture, etc. lt is important to 
identify the characteristic features for objects of 
interest in order to differentiate objects of different 
land cover types. Usually, it is based on a time 
consuming trial and error process. According to the 
selected training objects and features, 

SEparability and Threshold (SEaTH) algorithm 
calculates the proper features and thresholds that are 
best fit for the separation of interested classes. Based 
on astatistical approach, it calculates the features and 
thresholds for any number of given features and any 
number of interested classes, so it avoided the time- 
consuming trial and error process. This paper 
investigates the performance of object based image 
analysis for mapping land-covers in a mountainous 
forested area, using the features and thresholds 
identified by the SEaTH algorithm. 


2. METHODS 
2.1. Object based image analysis in eCognition 


Object based image analysis in eCognition (Definiens 
2006) comprises two parts: multi -resolution 
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segmentation and context-based classification. Multi - 
resolution segmentation allows generating image 
objects on an arbitrary number of scales considering 
homogeneity in colour and shape. 

Additionally, the created segments/objects are 
embedded into a hierarchical network in which each 
object knows its neighbouring objects in horizontal 
and vertical direction (figure 1) (Baatz and Schape 
1999). The classification is based either on a fuzzy 
nearest neighbour classifier or membership function. 


2.1.1. Multi-resolution image segmentation 


Multi -resolution image segmentation extracts objects 
in a number of hierarchical segmentation levels, and 
each subsequent level yields image objects of a larger 
average size by combining objects from a level below, 
which represents image information on different 
scales simultaneously (figure? ). 

In image segmentation, scale parameter determines 
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the maximum allowed heterogeneity for the resulting 
image objects. With a certain scale parameter, objects 
are formed based on “colour” and “shape” 
characteristics, and the relation between these two 
criteria is colour + shape = 1. “Colour” defines the 
weight the spectral values of the image layers 
contribute to image segmentation. Shape criterion is 
represented by smoothness and compactness. 
Smoothness is to optimize image objects with regard 
to smooth borders, and compactness with regard to 
compact borders (Baatz et al. 2004). 

Besides, for a given set of segmentation parameters, 
heterogeneous image data result in smaller image 
objects than homogeneous image data. 


2.1.2. The fuzzy classifier and hierarchical 
network classifier 


The classification uses either nearest neighbour 
classifier (NN) or membership function, or a means 


Image 
e object 
hierarchy 
nel level 


Figure 1: Image object hierarchy. Source: Definiens professional 5 user guide 


combination of both. Object features can be used 
either to describe membership functions, or to 
determine the feature space for NN (Baatz and 
Schape 1999). NN classifier is based on fuzzy logic, 
which replaces the two Boolean logical statements 
“true” and “false” by the continuous range of (0...1) 
where O means false and 1 means true and all values 
between 0 and 1 representa transition in between. 

To guarantee a definite hierarchy the segmentation 
procedures follow two rules: 1) object borders must 
follow borders of objects on the next lower level; 2) 
segmentation is constrained by the border of the 
object on the next upper level. Different hierarchical 
levels can be segmented using different data, an 
upper level, for example, can be built using thematic 
information, whereas a lower layer is segmented 
using remote sensing data. Classifying the upper level, 
each object can be analyzed based on the 
composition of its classified sub-objects. By means of 
this technique different data types can be analyzed in 
relation to each other. 


2.2. SEaTH algorithm 
SEaTH algorithm was proposed by Marpu et al. 


(2006) and Nussbaum et al. (2006), the algorithm 
identifies the characteristic features with a statistical 


approach based on training objects. lt measures the 
statistical distance (Jeffries-Matusita distance) to 
compare the feature separability of two interested 
classes. The measure result is in interval of [0, 2], in 
which O implies that the class distributions in these two 
features are completely correlated and 2 implies the 
complete uncorrelation. 

The features in which the class has the high 
separability (high J value) are the optimum features. 
To apply this method in the classification in eCognition, 
first of all, samples of the interested classes are 
selected and the statistics (feature values) of those 
samples are exported. 

The separability of the classes represented by J 
distance is then calculated, and the features with the 
high J distance are favorable features to differentiate 
between land-cover classes. The selected favorable 
features can then be used to build the rule base. Itis a 
parametric approach and the assumption is that all the 
random variables encountered during the processing 
can be approximated as Gaussian distributions. 


3. THESTUDY AREA AND DATA 
The study area is located in Michoacán state, central 


west of Mexico, covering an area of approximately 
58*60 km 2, within the longitude of 102” 00"W and 
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102” 32" W, and latitude of 19” 02” N and 19* 36" N 
(figure 2). 


Michoscán Sine 
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Figure 2: The study area. Left side of the figure are two sketch maps indicating Mexico and Michoacán state 
where the study area is located; right side is the false colour composite of Landsat-7 ETM+ image with red, 
green, and blue bands of 4, 5, and 7. 


The available data are comprised of a Landsat-7 
ETM+ image obtained on 16 February 2003, with 6 
bands and a spatial resolution of 30m, a mosaic of 25 
ortho-corrected photos taken in 1995 with 2 meters 
spatial resolution, and a land-cover map from the 
Forest Inventory generated in the year 2000. Eight 
land-cover types need to be classified: “orchards”, 
temperate forest', 'irrigated agriculture”, “rain-fed 
agriculture”, lava flow”, “tropical dry forest', “grass 
land”, and 'human settlement. 


4. RESULTS AND DISCUSSION 


Principal component analysis was carried out to the 
six Landsat-7 ETM+ spectral bands and three PCA 
bands were generated. NDV| was calculated based 
on bands 3 (red), and 4 (near infrared) and its value 
was calibrated into the range of 0-300. The three PCA 
bands, the calibrated NDVI, and the six Landsat 
spectral bands were used in the object based image 
analysis. Two segmentation levels were created by 
multi -resolution segmentation. In the first level image 
was segmented using both image layers and thematic 
layer which contains 'human settlement' information. 
Including the information of a thematic layer in the 
segmentation influences the generation of image 
objects. Thematic layers contain discrete information, 
which means that related layer values can carry 


different information that is defined in the attached 
attribute list. To be able to clearly define the affiliation 
of an object to a thematic class given by the themati c 
layer, it is not possible to create image objects which 
belong to different thematic classes. To ensure this, 
the borders separating different thematic classes are 
restrictive for further segmentation whenever a 
thematic layer is used. For this reason, thematic 
layers can not be given different weights but can 
merely be selected for use or not. In this case, of 
course, the thematic layer which contains the 
settlements information was given the weight one to 
include it in the segmentation. Meanwhile, the image 
layers can also be given to different weights. If switch 
the weights for all image layers to O, the produced 
image objects are based exclusively on thematic layer 
information and the obtained objects do not represent 
the averaging of layer values. The purpose of the 
segmentation in this level was to create only 
settlement and non-settlement objects. In order to 
group all the objects which are not the settlement 
objects into one object for facilitating the object 
allocation, an extremely large scale factor 10000 was 
used and the segmentation created “human 
settlement' objects and one non-settlement object 
(figure 3). This level acted as super level in class 
hierarchy. 
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Figure 3: Segmentation using the thematic layer 
(level 1) 


A second segmentation level was created with the 
three PCA bands and the calibrated NDVI with 
parameters of scale factor 25, colour factor 0.7, 
compactness 0.5. The selection of appropriate 
segmentation factors was based on trial and error with 
the segmentation procedure until a satisfactory 
pattern was found. In this case, the produced 
segments are smaller than the objects of interest. The 
importance of visual inspection in segmentation result 
evaluation was addressed in Pal and Pal (1993) and 
Benz (2004) and it is still a common practice in 
segmentation parameters selection. The settings of 
all segmentation parameters are very much 
dependent on the image data, study area, and the 
desired land cover types, which gives the 
segmentation black box characteristics, though 
several groups are working on approaches to 
overcome this limitation (e.g. Espindola et al. 2006). 

These two levels were differentiated using 
hierarchical membership function. The second level is 
the actually classification level (figure 4). The smaller 
objects (in level 2) are nested within larger 5 ones (in 
level 1) and inherit any characteristics of a larger 
image object. First the super-level was classified. Two 
child classes were created in this level: 'human 
settlement' and 'non human settlement”. Eight classes 
were created in the classification level: “human 
settiement”, “irrigated agriculture”, “rain-fed 
agriculture”, “grassland”, “tropical dry forest”, 
'orchards”, temperate forest”, and lava flow". "Human 
settliement' at this level was classified using the 
classification result from its super level “class related 
features, relation to super class 'human settlement . 
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The rest seven class use also class related features 
from the super level class 'non-settlement”: “class 
related features, relation to super class 'non- 
settlement”. Besides this class related feature, other 
features used were calculated using the SEaTH 
algorithm. For SEaTH analysis about 2-5% sample 
objects of each land cover class were selected as 
training data. Then 104 characteristic features 
including spectral related, shape, and texture features 
were selected for the training objects. Based on the 
training objects and selected features, SEaTH tool 
calculated, for every object class combination, the 
separability and thresholds for each of the 104 
features. Based on the separability value, the “best” 
among 104 features are selected for each object class 
combination, and the features with maximum 
separability are used for classification. The selected 
features and thresholds based on SEaTH results for 
the land-cover classes were shown in the figure 5. 

The calculated features and thresholds are presented 
in the “class description” in eCognition. A class 
description is composed of expressions which 
describe a class through fuzzy or crisp terms and are 
connected by logical operators. Logical operators are 
needed because usually there are multidimensional 
dependencies in the feature space and a logical 
combination of features is used to represent this 
condition. The most common operators are “and 
(min) and “or (max)” which are used in the following 
way: 





Figure 4: Segmentation by scale factor 20, 
colour 0.7, and compactness 0.5 (level 2) 
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Figure 5: The calculated features and thresholds for the land cover classification by SEaTH algorithm 


In the classification scheme, class 'irrigated 
agriculture” has three sub-classes and 'temperate 
forest' has five sub-classes. Using SEaTH algorithm, 
the training objects were selected for each of the sub- 
classes and SEaTH calculated the features and 
threshold for the sub-classes. In eCognition, after the 
classification of the sub-classes, they can be grouped 
together as the classes of interest. 

Most of the object primitives were classified correctly 
into the classes of interest. This classification does not 
produce fuzzy results. Due to the spectral variety of 
objects of interest, it is possible that the selected 
training samples do not always fully represent the 
characteristics of those objects, and thus there are 
objects that were not classified. If there are a large 
amount of objects unclassified, the process of 
selecting the training objects and calculating the 
features and thresholds should be re-done; if the 
unclassified objects are relatively few , the calculated 





Figure 6: Classification result of level 2. 


thresholds can be adjusted in order that the feature 
values of those objects are being included and thus 
they can be classified. 

First it needs to be clear which classes those 
unclassified objects should belong to, and then the 
object feature values are compared with those 
desired classes. However, adjusting the calculated 
threshold could cause the mis-classification of other 
objects and it should be done carefully. When very few 
objects remain unclassified, manual classification can 
be used to allocate them to the correct classes. In this 
case, there were only a reduced number of objects 
unclassified. By observing the characteristics of those 
objects and comparing them with the features and 
thresholds of the desired classes, the thresholds of 
the classes were adjusted and the objects were 
reclassified. 

The classification resultis shown in figure 6 and 7. 
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Figure 7: Legend of SEaTH classification result 
(level 2) 
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Accuracy assessment was carried out based on 600 
stratified random points. 75 random points were put to 
each of the eight classes of interest in a land use map 
(2000). These points were interpreted based on the 
ortho-corrected photographs, land use map 2000, 
data from field survey, and Landsat-7 ETM+ images. 
The accuracy assessment result was then put into an 
error matrix based on the ground data and 
classification data. Indices such as the overall 
accuracy, the producer's and user's accuracies were 
calculated and it obtained an overall accuracy 78.7%. 


5. CONCLUSIONS 


Image objects are characterized by a variety of 
spectral, spatial, texture, and contextual features. 

The usually trial and error process of searching the 
proper features is time consuming and thus hinders 
the utilization of the strength of object based image 
analysis. Separability and Threshold (SEaTH) 
algorithm is able to evaluate statistically any number 
of given features for any number of object classes of 
interest and to identify the optimum features and 
thresholds for the separation of objects of interest. 
The working process of SEaTH algoritam can be 
generally described as the following: first, training 
objects for the classes of interest were selected. lf a 
class has several subclasses, training objects should 
be selected for each of the sub-class. Then the feature 
values of those training objects are exported and the 
desirable features and thresholds are calculated 
based on the separability. For each class, SEaTH 
compares itother classes using each ofthe selected 
features. The features and thresholds were listed in 
the sequence of the separability values. Based on the 
calculation, features and thresholds with the highest 
separability-the highest Jeffries-Matusita (J) value, 
which is in the range of (0, 2), are selected for the 
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RESUMEN 


La llanura aluvial actual del Paraná, conforma una 
llanura viva, compleja, dendrítica, la porción entre 
Santa Fe y Rosario con una compleja red de drenaje 
anastomosado con una alta densidad de cuerpos de 
agua grandes, irregulares e interconectados, con 
islas de mayor tamaño hacia el sur, y con áreas 
geomorfológicamente más activas son las ubi cadas 
en Oo cerca de los canales. El área de estudio, 
delimitada desde la diagonal Paraná/Santa Fe 
(31%40"S) hasta Puerto Gaboto/Arroyo Las Ceibas 
(32305), corresponde a un área de ambientes 
lóticos, lénticos, y de islas con diferente tamaños y 
dinámicas de inundación. El sistema hidrológico del 
Río Paraná en este tramo ejerce a escala regional un 
efecto de control. A escala local los principales 
procesos de transformación están ligados a las 
historias de uso del suelo. La Teleobservación se 
presenta como una herramienta de análisis 
apropiada, y una fuente de datos adecuada para el 
estudio de las características y dinámicas 
morfológicas y evolutivas de la cobertura del suelo de 
ambientes de la llanura aluvial del Paraná. Se ha 
definido como objetivo general del trabajo ensayar 
técnicas de procesamiento de imágenes para la 
caracterización de la cobertura de suelo en la llanura 
aluvial del Río Paraná de acuerdo a su geomorfología, 
topografía e hidrología. Se trabajó sobre una ventana 
de la imagen Landsat TM5 del 13 de febrero de 2007, 
Path/Row 227/082. Se utilizaron también cartas 
topográficas del Instituto Geográfico Militar (IG M) 
escalas 1:100.000 y 1:50.000. Se generó cartografía 
temática correspondiente a geomorfología, 
topografía e hidrología. Los resultados obtenidos 
indican que el área estudiada presenta una alta 
heterogeneidad espacial (y temporal) dada 
fundamentalmente por la dinámica de los procesos 
geomorfológicos y fluviales, modeladores principales 
del paisaje. Se diferenciaron grandes unidades 
ambientales asociadas a diferentes procesos 
modeladores, con caracte-rísticas y evolución 
propias. Se reconocen sistemas de islas 
evolutivamente clasificadas en complejas y simples, 
correspondientes a depósitos de cauce, con 


sedimentación de mayor granulometría y vinculadas 
directamente al régimen hidrológico de alta energía. 
Se identificaron así mismo, sistemas más estables, 
de tierras altas (porción mas externa de la Dorsal del 
Paraná) y bajas (sistema Coronda y parte sur de la 
porción superior del Complejo Deltaico). 

Palabras clave: unidades de ambiente, Río Paran á, 
coberturas, geomorfología, vegetación, hidrología. 


ABSTRACT 


The current floodplain of Paraná, forms a living, 
complex and dendritic plain, the portion between 
Santa Fe and Rosario with a complex network of 
drainage anastomosed with a high density of large 
irregular and interconnected bodies of water, with 
larger islands towards the south, and with areas 
geomorphology more active are located on or near the 
canals. The area of study, bounded from the diagonal 
Parana / Santa Fe (31 * 40'S) to Puerto Gaboto / The 
Arroyo Ceiba (32 * 30'S), corresponds to an area of 
lotic and lénticos environments, and islands with 
different sizes and dynamics of flooding. The 
hydrological system in Rio Parana in this stretch 
exercises at regional level an effect of control. At local 
level major transformation processes are linked to the 
stories of land use. Remote sensing ¡is presented as 
an appropriate analysis tool, and a source of data 
suitable for studying the characteristics and 
morphological and evolutional dynamics of ground 
cover in environments of the floodplain from Parana. lt 
has been defined as a general objective in this study 
to try out image processing techniques for the 
characterization of covering ground on the floodplain 
of the Parana River according to its geomorphology, 
topography and hydrology. We worked on a window of 
the Landsat image TM5 from February 13, 2007, Path 
/ Row 227/082. It was also used topographic charts 
from the military geographic institute (IGM) scales 
1:100.000 and 1:50.000. lt was generated thematic 
cartography for geomorphology, topography and 
hydrology. 

The results obtained indicate that the studied area has 
a high spatial heterogeneity (and temporary) given 
mainly by the dynamics ofthe fluvial and 
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geomorphologic processes, main modell ers of the 
landscape. Several environmental units associated 
with different processes modellers with own 
characteristics and evolution were differed. lt can be 
recognized systems of islands evolutionarily 
classified in complex and simple, for deposits channel, 
with sedimentation of greater particle size and link 
directly to the high energy hydrological regime. 
Moreover, more stable systems were identified, 
higher land (largest portion of the Dorsal Parana) and 
low land (Coronda system and southern part of the 
upper portion of the deltaic Complex). 

Key words: environmental units, Paraná River, 
coverage, geomorphology, vegetation, hydrology. 


INTRODUCCIÓN 


Los sistemas de humedales están estrechamente 
vinculados al régimen hidrológico (Neiff, 1996) 
dentro del paisaje de la llanura de inundación, la 
morfología y los depósitos sedimentarios están 
determinados por el régimen de flujos y el grado de 
conectividad entre sistemas lóticos y lénticos que lo 
constituyen (Drago, 1990; Drago et al., 2003; Drago, 
2007; Iriondo, 2007). Así, el flujo lateral de agua 
durante los períodos de inundación, se convierte en la 
fuerza principal que organiza y regula el 
funcionamiento de los humedales riparios (Malanson, 
1993). 

Diversos estudios realizados sobre la llanura aluvial 
del Río Paraná (Iriondo, 1972, 1991, 2004, 2007; 
Iriondo y Drago, 1972; Drago, 1973, 1981, 1990, 
2007; Iriondo y Paira , 2007) reconocen la 
heterogeneidad y complejidad regional, la dificultad 
de consensuar terminología dada las diferentes 
escalas a las que se expresan diversos procesos, y la 
estrecha vinculación de los procesos geomor- 
fológicos e hidrológicos , que convergen en la 
geomorfología fluvial. La conectividad entre estos 
subsistemas se ha utilizado para clasificar los 
ambientes acuáticos (Drago 1973, 1990; Paira y 
Drago, 2006). 

Oakley et al. (2005) proponen una regionalización del 
Corredor Fluvial Paraguay-Paraná, donde el sector 
caracterizado en presente trabajo es incluido en la 
“Gran Unidad Natural Paranaense”, y en la “Gran 
Unidad Natural Delta del Paraná”. 

Para el Valle de aluvial del Paraná, el primer modelo 
interpretativo sobre la vegetación utilizando métodos 
fitosociológicos fue realizado por Morello (1949) en 
las islas del río Paraná situadas al sur de nuestra área, 
frente a la cuidad de Rosario. Lewis y Franceschi 
(1979) y Franceschi y Lewis (1979) caracterizaron la 
vegetación y elaboraron un concepto dinámico de las 
mismas. Aceñolaza et al. (2005) y Aceñolaza y 
Marchetti (2005) determinaron unidades 
homogéneas de vegetación para islas del Bajo 
Paraná ; bosques de barranca, bosques simples de 
albardones marginales y mixtos de albardones 
internos, unidades de ambientes intermedios, 
vegetación lacunar y pastizales degradados. 

Para el extremo superior del Bajo Paraná, Neiff (1999, 
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2002) y Casco (2004) estudiaron la distribución de la 
vegetación en relación con la hidrología. Para el 
sector del Delta del Río Paraná, Malvárez (1999) 
integró unidades de paisaje, suelos y vegetación; 
correspondiendo al sur del sector estudiado en el 
presente trabajo la unidad de “bosques, praderas y 
lagunas de llanuras de meandros”. 

Los autores previamente citados coinciden en que las 
fisonomías arbóreas ocupan las porciones altas, las 
herbáceas porciones intermedias y las acuáticas las 
porciones más bajas del gradiente topográfico; 
estableciendo que gran parte de la vegetación fluvial 
se distribuye en un amplio rango de condiciones de 
hábitats, expresando así su plasticidad a las 
variaciones ambientales. 

En este contexto, se plantea como objetivo general 
del presente trabajo el ensayo de técnicas de 
procesamiento de imágenes para la caracterización 
de la cobertura de suelo en la llanura aluvial del Río 
Paraná de acuerdo a su geomorfología, topografía e 
hidrología. El desarrollo del mismo permitirá generar 
herramientas para estudios científicos y técnicos 
regionales, así como para la planificación e 
implementación de estrategias de gestión del 
territorio coherentes con sus características 
ambientales. 


EL ÁREA DE ESTUDIO 


El área de estudio queda delimitada a la región 
superior del complejo Deltaico del Río Paraná, desde 
la diagonal Paraná/Santa Fe (31%40"S) hasta Puerto 
Gaboto/Arroyo Las Ceibas (32305) (Figura 1). 

Las características climáticas del área la ubican 
dentro de un clima templado/cálido y húmedo. 

Las temperaturas medias anuales se encuentran 
cercanas a los 19%C. El régimen pluviométrico es de 
900/1000 mm, con precipitaciones que se registran 
principalmente en el período de octubre a abril (73%) 
(Rojas et. al, 1987). Dentro del contexto regional, la 
llanur a aluvial del Paraná contribuye a morigerar las 
condiciones climáticas del corredor, lo que estaría 
dado por la influencia de un gran volumen de agua 
fluyendo en dirección norte -sur (Iriondo, 2007). 
Desde un punto de vista geomorfológico el área de 
estudio forma parte de la Llanura aluvial actu al del 
Paraná (Iriondo, 1991), limitada hacia el Este por el 
bloque mesopotámico y por el Oeste por las regiones 
del Chaco y la Pampa (Iriondo y Paira, 2007) ; y se va 
ensanchando a medida que el río se extiende hacia el 
sur. A la altura de la Ciudad de Diamante el cauce se 
separa de la margen izquierda (Iriondo, 1972; Drago 
1973; Iriondo 2007). Conforma una llanura viva, 
compleja, dendrítica (Drago, 1973), con su parte alta 
en subsidencia, la porción entre Santa Fe y Rosario 
con una compleja red de drenaje anastomosado con 
una alta densidad de cuerpos de agua grandes, 
irregulares e interconectados (Paira y Drago, 2007) y 
con las islas de mayor tamaño hacia el sur (Iriondo, 
2007). Las áreas geomorfol ógicamente más activas 
son las ubicadas en o cerca de los canales (Paira y 
Drago, 2006). 
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Iriondo (1972), Iriondo y Drago (1972) y Drago (1973) 
identificaron depósitos producidos en el cauce y fuera 
del mismo por inundación (áreas más antiguas). Los 
primeros conforman llanuras de bancos (con bancos 
con sus ejes mayores paralelos a la dirección del flujo; 
Drago, 1990) , derrames y llanuras de meandros 
(Iriondo, 2007) . Las islas evolucionan de formas 
simples, hacia formas con lagunas en su interior, y 
bancos en áreas vecinas que posteriormente se 
adosan por acumulación de sedimento en las bocas 
de los cauces, resultando islas más complejas y con 
menor influencia directa del curso de agua adyacente 
(Drago, 1973). 

Los depósitos de Inundación, forman llanuras con 
avenamiento impedido y dendrítico (con lagunas 
circulares, de poca profundidad y menos extensión, 
cauces divagantes y alta densidad hidrográfica) 
bañados, pantanos y albardones (Drago, 1973). Otro 
tipo de depósitos (derrames) corresponde al delta del 
río Coronda, que sedimenta lateralmente dentro de la 
llanura aluvial, con llanuras de meandro y de 
avenamiento impedido vinculadas a este sistema 
(Iriondo, 1972, 1991). Como depósitos antiguos 
Iriondo (2007) reconoce áreas adyacentes elevadas 
formando terrazas, superficies con relieves bajos 
mayormente no visibles. Ubicadas en el sector 
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oriental de la llanura, y recorriendo 1200 km. y con 30 
a 50 mts. de elevación se encuentran las barrancas, 
nicho para la vegetación tropical en latitudes más 
templadas (Aceñolaza et al., 2005; Iriondo, 2007). 

El patrón hidrológico del río Paraná a esta latitud, 
posee un régimen pulsátil anual de inundación. La 
época de estiaje (bajante) se produce en agosto - 
septiembre, mientras que el pulso de inundación 
suele darse sobre fines del verano a principios del 
otoño (Drago, 1990; 2007). En esta región, el río 
Paraná contiene un mosaico heterogéneo de cuerpos 
de agua lóticos y lénti cos que cambian de forma y 
área de acuerdo a las dinámicas fluviales (Paira y 
Drago, 2007). El régimen hidrológico actúa como 
modelador formando canales entrelazados con 
ensanchamientos y estrangulamientos (Drago etal., 
2003), con procesos de erosión y sedimentación 
lateral que originan la formación y desaparición de 
islas (Drago, 1990; Iriondo 1991; Iriondo, 2007 ; 
Iriondo y Paira, 2007). En los períodos de inundación 
suele haber erosión en las partes anteriores de las 
islas y sedimentación en las posteriores, causando la 
variación de formas y migración de las islas aguas 
abajo (Drago, 1990). 
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Figura 1: Área de estudio. 
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MATERIALES Y MÉTODOS 


Para caracterizar el área de estudio, se realizó una 
revisión de los trabajos existentes, su cartografía e 
interpretaciones. Á partir de esos datos y del análisis 
de imágen es satelitales, se estructuró una nueva 
serie de mapas correspondientes a las unidades 
geomorfológicas, hidrología y de paisaje. 

Las interpretaciones se realizaron a partir del análisis 
visual y digital de imágenes satelitales. Para el 
procesamiento de imágenes satelitales se trabajó 
sobre una ventana de la imagen Landsat TM5 del 13 
de febrero de 2007, Path/Row 227/082. Se utilizaron 
también cartas topográficas del Instituto Geográfico 
Militar (IGM) escalas 1:100.000 y 1:50.000. La 
ventana de imagen fue pre - procesada, georrefe- 
renciándose a partir de puntos obtenidos sobre 
imágenes de alta resolución del Google Earth y de 
coordenadas propias obtenidas con GPS. El método 
utilizado fue el de vecino más próximo para preservar 
los valores originales de la ventana correspondiente a 
la imagen seleccionada de 2902 por 4516 píxeles. No 
se efectuó corrección radiométrica de los datos ya 
que el presente estudio no fue de carácter multi 
temporal. Para el proceso de clasificación fueron 
considerados los ND originales de la imagen. 

Para la identificación de unidades de paisaje y la 
vegetación característica de las mismas, se realizó 


una clasificación no supervisada (ISODATA) (Tou etal. 


1973; Chuvieco, 2002) con 15 clases. Posterior- 
mente se procedió a asignar cada una de las clases 
generadas a las diferentes unidades de vegetación y 
ambiente definidas en la leyenda. Para la correcta 
asignación se utilizó información propia de terreno e 
imágenes de alta resolución disponibles desde el 
Google Earth. Estas clases fueron luego reagrupadas 
en 4 clases en base a los datos relevados en te rreno. 
Para la determinación de cuerpos de agua y su 
posterior discriminación entre cuerpos con mayor 
presencia de sedimentos en suspensión 
(frecuentemente lóticos ) y de aquellos con menor 
proporción de sedimentos en suspensión (lénticos ), 
se utilizó un proceso de segmentación basado en los 
niveles digitales de los píxeles. Se determinaron 
sobre la imagen los umbrales mínimos y máximos 
para los diferentes cuerpos de agua y a partir de esos 
rangos de valores se procedió a segmentar la imagen 
en base a los niveles digitales de los canales 
correspondientes al azul, infrarrojo cercano e 
infrarrojo medio. 

El criterio de segmentación sobre los niveles digitales 
para la identificación de todos los cuerpos de agua 
sobre la imagen fueron: Azul (Banda 1 Landsat 5TM ) 
> 45, Infrarrojo Cercano (Banda 4 Landsat5TM ) < 32, 
Infrarrojo Medio (Banda 5 Landsat5TM ) < 22. 

Para la identificación de los cuerpos de agua Lénticos 
sobre la imagen se utilizaron los siguientes criterios 
de segmentación sobre los niveles digitales de la 
imagen: Azul (Banda 1 Landsat 5TM ) > 45 y < 65, 
Infrarrojo Cercano (Banda 4 Landsat 5TM ) < 20 , 
Infrarrojo Medio (Banda 5 Landsat5TM)< 15. 

Por último, se visitaron sitios representativos de las 
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diferentes unidades de vegetación y ambiente en el 
área de estudio, para los cuales se tomaron 
composiciones florísticas , cobertura y tipo de estrato 
vegetal presente. Para cada uno de los sitios 
relevados a campo se registraron coordenadas GPS 
y setomaron fotografías. 


RESULTADOS Y DISCUSIÓN 


A partir del análisis de la cartografía vigente, y del 
propio análisis, se presentan una serie de mapas 
correspondientes a las interpretaciones de las 
unidades Geomorfológicas , hidrología y de paisaje. 
Iriondo (2007) distingue 3 tipos de procesos de 
sedimentación (Figura 2): correspondientes a 
depósitos de cauce (1), antiguos (2), y de inundación 
(3). Los depósitos de cauce se distribuyen cercanos a 
los cursos de agua principales , corresponden a 
sedimentación de textura más gruesa generada por la 
mayor energía del río y su corriente turbulenta. Son 
depósitos dinámicos con continua remodelación por 
erosión y deposición de nuevos sedimentos 
transportados. Los depósitos de inundación siguen 
una línea paralela al margen por el cual el cauce se 
desborda, reconociéndose depósitos de inundación 
asociados al sistema río Paraná (Fig.2) y otros al río 
Coronda (Fig. 2). 

Estos depósitos de textura más fina, sedimentan 
durante o después del pico de inundación, por 
disminución de la corriente del agua y fijación por la 
vegetación. 





Figura 2: Mapa de depósitos elaborado a 
partir de Iriondo (1972). 
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Figura 3: Cartografía con visualización de 
sistemas. 


En la figura 3 se identifican 5 sistemas que confluyen 
en el área a partir de la interpretación de las 
cartografías existentes y las propias. Hacia el centro 
del corredor fluvial se diferencia un sector de islas del 
cauce principal e islas adyacentes (Fig. 3). Las islas 
adyacentes al cauce son depósitos que, en algunos 
casos, conforman la desembocadura de ríos y 
arroyos de la Provincia de Entre Ríos sobre el Río 
Paraná. Algunas de estas islas fueron descriptas por 
Iriondo (2007) como “depósitos antiguos”. En el 
sector central del área de estudio se encuentra una 
unidad que corresponde a lo que hemos denominado 
dorsal del Paraná, que representa una zona de mayor 
elevación topográfica que recorre el margen 
occidental del río, con un gradiente de mayor a menor 
hacia el oeste, estructura similar a la observada en 
latitudes menores (sistema Saladillo). El sector 
occidental corresponde a la cuenca del río Coronda 
(Sistema Coronda), que representa una estructura 
con derrames secundarios laterales. Hacia el sur se 
ha identificado un sector construido a ambas 
márgenes del actual cauce principal del Paraná, que 
se extiende desde la altura de la localidad de General 
Alvear (Entre Ríos), formando un abanico que va 
ensanchándose hasta alcanzar más de 30 km a la 
altura de la desembocadura del Río Coronda al 
Paraná y constituye la porción superior del complejo 
Deltaico (Fig. 3, 5). 

La figura 4 revisa los criterios previamente 
mencionados , identificando 4 categorías. En primer 
lugar se diferencian tierras altas y bajas , definidas de 
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acuerdo al gradiente topográfico. Las tierras altas 
dentro de la llanura de inundación, han sido 
construidas mayoritariamente por el sistema del 
Paraná, y su hidrología está determinada por el 
patrón de pulsos de este sistema fluvial. Por otro lado 
las tierras bajas responden a dos dinámicas fluviales: 
la del río Coronda en el segmento occidental del área 
de estudio y la del Paraná en el oriental. Las otras 
categorías corresponden a las islas, las que 
presentan diferentes grados de evolución, 
diferenciándose en complejas (con cuerpos de agua 
en su interior, con menor conectividad al cauce 
principal) e islas simples. Para el área considerada, 
los cuerpos de agua (Figura 5) suman 111.200 has de 
las que 66.900 has. corresponden a sistemas con 
mayor contenido de sedimentos en suspensión. Esto 
se asocia a una mayor movilidad de agua ya sea por 
tratarse de agua corriente (sistemas lóticos) y/o 
grandes lagunas con circulación asociada a remoción 
por vientos (Coronda, Campo Grande, etc.). De las 
44.300 has de cuerpos de agua con baja cantidad de 
sedimentos en suspensión, 18.500 has corresponden 
a la Provincia de Entre Ríos, aunque quedan 
excluidas grandes áreas de lagunas de la provincia 
de Santa Fe por haber sido categorizadas como “con 
sedimentos en suspensión”. En líneas generales 
puede observarse la presencia de grandes lagunas 
asociadas al margen derecho y lagunas de tamaños 
menores sobre el margen izquierdo del área de 
trabajo. La relación de superficies entre laguna/tierra, 
es mayor sobre el margen santafecino que en la 
porción que se encuentra al este del cauce del río 
Paraná. 














Figura 4: Caracterización de áreas por alturas relativas e 


islas por grado de desarrollo. 


2 VOL. 26 N” 1 


JUNIO 2008 TELEOBSERVACIÓN APLICADA A LA CARACTERIZACIÓN... 














Figura 5: Cuerpos de agua. 1: Con mayor proporción de 
sedimentos en suspensión; 2: con menor proporción de 


sedimentos en suspensión. 


En relación al análisis de unidades de vegetación, el 
análisis espectral de las imágenes satelitales, 
muestra un patrón general que presentan fisonomías 
arbóreas dominando las áreas altas, fisonomías 
herbáceas altas o arbustivas en topografías 
intermedias y herbáceas bajas en las áreas planas 
más bajas. En este sentido, la vegetación del área 
puede agruparse en cuatro clases: bosques; 
arbustiva/herbácea alta; herbácea baja y acuática. 

La vegetación de bosque esta conformada por cuatro 
tipos diferentes: bosques de barranca, bosques 
simples de albardones marginales, bosques simples 
de albardones internos y bosques mixtos de 
albardones internos. Asimismo, es frecuente 
encontrar áreas cubiertas por bosques de baja 
densidad (tipo sabana ) en áreas ecotonales entre 
estos bosques y los pastizales/arbustales vecinos . 
Un patrón común para esta fisonomía, es una 
importante relación largo/ancho, condicionada por las 
características morfológicas de los sitios de 
colonización, generados por la actividad fluvial. 

Los bosques de barranca recorren el margen 
izquierdo del área de estudio, cubriendo las 
barrancas entrerrianas del río Para ná y no sufren 
inundaciones . En los mismos, se conjugan tanto 
especies de la llanura de inundación propiamente 
dicha como especies del Espinal, confiriéndole 
especial singularidad algunas pocas especies 
subtropicales que llegan a estas latitudes gracias a su 
sola presencia y que no forman parte de las dos 
unidades anteriormente mencionadas para esta 
latitud (Aceñolaza et al. 2005) . En el área de estudio 
esta unidad presenta variabilidad florística asociada a 
las características geomorfológicas de la barranc a 
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(pendiente, exposición, forma) y a sus respectivos 
estadíos sucesionales . Entre los primeros se destaca 
la presencia de Prosopis nigra, Capparis tweediana, 
Celtis spp. Porlieria microphylla, Dodon aea viscosa, 
etc. En los bosques de mayor cobertura, y en una 
etapa sucesional posterior, aparecen enriqueciendo: 
Phytolacca dioica , Myrcianthes cisplatensis, 
Ruprechtia laxiflora, Myrcianthes pungens, Myrsine 
laetevirens, Fagara hyemalis , Coccoloba 
argentinensis. El estrato herbáceo es muy rico en 
especies con numerosas especies de gramíneas y 
latifoliadas y con dominancia de unas sobre otras 
dependiendo del grado de cobertura vegetal. 

Sobre el valle aluvial propiamente dicho, 
encontramos los bosques simples de albardones 
marginales, los que están constituidos 
mayoritariamente por stands monoespecíficos y 
coetáneos de Salix humboldtiana y se distribuyen a lo 
largo de los albardones marginales de los ríos y 
arroyos principales con una baja recurrencia de 
inundaciones. También se encuentran bosques con 
similares características de Tessaria integrifolia, o 
aquellos constituidos conjuntamente por estas dos 
especies. Sobre los albardones internos de las islas 
pueden encontrarse dos tipos de bosques, los 
bosques mixtos, ricos en especies y asociados a 
albardones más altos (marginales de arroyos de alta 
energia) y los bosques simples asociados a áreas 
más planas o albardones más bajos. Los primeros, 
son las comunidades de mayor riqueza específica de 
las islas, y en él se conjugan algunas especies de los 
bosques simples de albardones marginales 
acompañadas con Nectandra angustifolia, Albizia 
inundata, Inga verna , Myrsine laetevirens y más de 
40 especies del sotobosque. Los segundos, son 
coetáneos y monoespecíficos de Sapium 
haematospermum o Albizia inundata. Estos bosques 
poseen una riqueza menor y proporcionalmente un 
mayor número de especies hidrófilas. 

La vegetación Arbustiva/herbácea alta conforma una 
unidad caracterizada por una serie de comunidades 
vegetales que se distribuyen por lo general en la 
porción alta de las medias lomas limitando con las 
unidades boscosas de isla y alternando en algunos 
casos con bosques abiertos de Salix humboldtiana o 
Acacia caven con fisonomía de sabana. Suelen tener 
suelos con mayor contenido de limo y arcilla que 
aquellos de los albardones marginales. Como 
comunidad herbácea alta, el pajonal de Panicum 
prionitis es una de las más representativas de la isla 
por su extensión espacial; asociada al disturbio del 
fuego e inundación recurrente que le permite 
mantener su estructura y composición . La vegetación 
Herbácea baja corresponde a una unidad que surge 
en época de estiaje sobre las áreas cedidas por las 
lagunas. Esta área posee un alto dinamismo de 
cambio intra-anual ya que puede encontrarse en 
situaciones alternadas de inundación y tierra firme 
que no permite la instalación del pajonal ni la 
subsistencia de plantas acuáticas sumergidas o 
flotantes. En esta circunstancia, estas amplias áreas 
son cubiertas por Alternanthera philoxeroides, 
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Ludwigia peploides, Paspalum spp., Polygonum spp. 
e Hydrocotyle bonariensis entre otras. Esta unidad 
posee alta productividad, siendo un área sometida a 
intenso pastoreo y quema por lo que se pueden 
encontrar situaciones diversas en cuanto a su 
composición florística, estado sucesional y biomasa. 
En este ambiente, se encontró que las comunidades 
herbáceas se regeneran rápidamente luego de 
disturbios catastróficos , como los de inundaciones 
extraordinarias, por lo que en su conjunto son 
resilientes a este tipo de eventos (Franceschi et al., 
2000). 

Por último, la vegetación Acuática se caracteriza por 
distintas comunidades flotantes , arraigadas y 
comunidades de plantas sumergidas. Prácticamente 
toda el área está permanentemente inundada, o se 
inunda con alta periodicidad, por lo que estas 
comunidades pueden cambiar en su importancia 
relativa por el cambio de extensión de las superficies 
que abarcan, siendo mayores en las zonas internas 
y/o bajas de las islas, y menores en los márgenes de 
los cursos de agua de alta energía (albardones). 
Entre la vegetación acuática arraigada, encontramos 
Nymphoides indica, Victoria cruziana, Polygonum 
acuminatum, Sagittaria montevidensis, Ludwigia 
peploides y Enhydra anagallis entre otras. Las 
comunidades acuáticas flotantes están constituidas 
por Salvinia biloba, Pistia stratiotes, Limnobium 
spongia, Azolla filiculoides, Wolffia spp., Lilaeopsis sp. 
Eichhornia azurea y E. crassipes. Por ultimo, 
encontramos una serie de especies vasculares que 
conforman comunidades acuáticas sumergidas. 
Estas se encuentran en madrejones y lagunas 
Ceratophyllum demersum, Utricularia gibba, Elodea 
callitrichoides, Egeria spp., que constituyen tapices 
subacuáticos en ambientes con bajo contenido de 
sólidos en suspensión como madrejones, lagunas y 
arroyos de baja energia. Integrando los resultados 
hasta aquí obtenidos se pueden plantear los 
siguientes puntos: 

El área estudiada presenta una alta heterogeneidad 
espacial (y temporal) dada fundamentalmente por las 
dinámicas de los procesos geomorfológicos y 
fluviales, modeladores principales del paisaje. Esta 
heterogeneidad puede ser interpretada como una alta 
capacidad de resiliencia de los sistemas que 
conforman las unidades ambientales. 

Los flujos de materia, energía y especies indican las 
direcciones de estos procesos modeladores. Así, la 
formación y desaparición de sedimentos en islas del 
cauce principal responde a la dirección del río Paraná 
(“migración de islas” en sentido N-S) mientras que en 
otros sectores el sentido de la deposición de 
sedimentos es perpendicular al cauce principal (flujo 
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lateral: Ej. Sistema Coronda). La alteración de estos 
flujos (y por lo tanto de la conectividad) significa la 
transformación de estos sistemas ecológicos. 

Se reconocen sistemas de islas, complejas (en su 
mayoría) y simples, correspondientes a depósitos de 
cauce, con textura de mayor granulometría y 
dinámica mas vinculada al régimen hidrológico. 
También se reconocen los sistemas más estables, de 
tierras altas (porción mas externa de la Dorsal del 
Paraná) y bajas (sistema Coronda y parte sur de la 
porción superior del Complejo Deltaico). 

La presencia del Sistema Coronda y la Dorsal del 
Paraná, genera una mayor proporción de cuerpos de 
agua sobre la margen occidental del río . La 
vegetación presenta patrones definidos (aunque 
acotados espacialmente) de distribución, asociados a 
gradientes topográficos (e indirectamente de 
conectividad fluvial). 

Por último, a los fines de la planificación y ordenación 
del territorio estudiado será de importancia tener en 
cuenta las siguientes consideraciones: 

Se diferencian grandes unidades ambientales 
asociadas a diferentes procesos modeladores, con 
características y evolución propias. Es importante 
reconocer las dinámicas (e interconexión) de cada 
unidad. 

No se observan sistemas altamente antropizados, o 
con fuerte alteración de sus elementos y procesos, 
por lo que puede considerarse un área apropiada 
para concentrar esfuerzos de conservación, 
aumentando así la representatividad en el Sistema de 
Áreas protegidas. El área estudiada está 
comprendida por las jurisdicciones de la Provincia de 
Entre Ríos (margen 

oriental) y la de Santa Fe (margen occidental). Estos 
límites políticos no deberían ser los únicos 
considerados a los fines de la gestión, ya que las 
unidades identificadas no responden a los mismos . 
Es importante en este punto reconocer la necesidad 
de consenso de las políticas territoriales de las 
diferentes administraciones. 
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RESUMO 


Este estudo tem como objetivo o desenvolvimento de 
uma metodologia para a classificacáo automática de 
manchas urbanas contínuas e dispersas, em regides 
com altos índices de urbanizacáo no Estado de Sáo 
Paulo. Busca-se desenvolver uma metodologia que 
permita repetir o procedimento em outras regides do 
país, com o aumento da exatidao do mapeamento, 
assim como a reducáo do tempo e a subjetividade. 
Para este fim utilizou-se a classificacáo Orientada a 
Objeto em imagens LANDSAT relativas a 1985 e 
1999. Este procedimento consiste na segmentacáo 
multiresolucáo das imagens e na classificacáo 
baseada na lógica fuzzy. Os resultados mostraram 
que a classificacáo orientada a objetos tem um bom 
desempenho para mapeamento de áreas urbanas 
com um índice kappa de 0,70 para o Município e 0,55 
para a área de uso dominantemente urbano. 


Palavras-chave: Análise orientada a objeto, 
urbanizacáo dispersa, sensoriamento remoto. 


ABSTRACT 


This study attempts to develop an automatic 
approach to map urban areas and urban sprawl in 
high urbanization regions of Sáo Paulo State. The 
objective is to achieve a methodology for the increase 
of precision and reduction of processing time, to be 
applied in other regions of Brazil. Object Oriented 
Image Analysis with eCognition software and 
LANDSAT images of 1985 and 1999 were used. The 
object-oriented approach is an areal analysis method 
based essentially on segmentation and classification 
outputs. Primitive objects are created on given space 
scale levels with different resolution, and 
classification rules are applied based on a fuzzy rule 
decision tree classifier. The results show that the 
Object Oriented Image Analysis have a good 
performance to detect urban features with kappa 
index of 0.70 (county) and 0.55 (principal urban area 
of county). 


Palavras-chave: object oriented image analysis, 


urban sprawl, remote sensing, 
INTRODUCÁO 


As cidades brasileiras apresentam crescimento 
rápido dado ás elevadas taxas de crescimento natural 
da populacáo e migracáo rural urbana, quando 
comparadas com países com alto desenvolvimento 
tecnológico. O monitoramento deste crescimento 
torna-se difícil, e oneroso, se realizado com técnicas 
convencionais de obtencáo de dados como 
levantamento de campo, levantamentos aero- 
fotogramétricos tradicionais embora este último 
tenha tentado ajustar-se paratornar-se competitivo. 
Desde o lancamento de satélites Landsat específicos 
para levantamento de recursos naturais iniciado em 
1972 tem-se através de suas imagens nova 
ferramenta para realizacáo de estudos urbanos 
devido a repetitividade dos dados e visáo sinótica 
fornecida pelas imagens associada á extensáo da 
área abrangida por uma única cena. 

Em um primeiro momento o uso de imagens de 
satélite direcionou-se para estudos de crescimento 
urbano horizontal com o mapeamento de mancha 
urbana, e estes apresentavam precisáo compatível 
com estudos regionais em funcáo da resolucáo dos 
dados e número de bandas espectrais dos sensores. 
As áreas de solo exposto preparadas para agricultura 
ou áreas de mineracáo, localizadas na franja rural- 
urbana, eram em sua quase totalidade classificadas 
como urbana. 

Com o avanco tecnológico no desenvolvimento de 
novos sensores resultando na melhoria da resolucáo 
espacial, radiométrica e espectral dos dados surgem 
novas possibilidades de uso de imagens orbitais para 
estudos urbanos mais detalhados. Ocorre também 
um avanco na melhoria da precisáo das informacdes 
geradas. 

Devido ao grande número de satélites existentes, 
aqueles desenvolvidos por órgáos governamentais e 
os comerciais, gerando grande oferta de dados, os 
custos na aquisicáo dos dados para estudos 
apresentam-se em declínio facilitando o acesso a 
estes dados. 

O lancamento dos satélites da série CBERS 
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resultante da cooperacáo China - Brasil abriu novas 
perspectivas para uso sistemático das imagens 
orbitais em estudos urbanos, principalmente em um 
país de dimens0es continentais como o Brasil. A 
melhoria na resolucáo espacial prevista para os 
novos satélites é importante para monitoramento das 
áreas urbanizadas e estudos do uso do solo urbano 
que além da complexidade os alvos urbanos tém 
pequena dimensáo associada a poucos metros. 

As informacóes retiradas das imagens a um baixo 
custo e em curto espaco de tempo poderáo ser 
integradas as informacóes obtidas por meio dos 
censos demográficos, pesquisas Institucionais 
realizadas, dados levantados por prefeituras, a partir 
do uso de Sistemas Geográficos de Informacáo e 
utilizados náo somente para o planejamento racional 
do uso de solo como também para a verificacáo da 
implantacáao das medidas estabelecidas pelos 
estudos desenvolvidos para elaboracáo de Planos 
Diretores. Alei, portanto, passará da teoria a prática. 
Com Oo avanco das técnicas de interpretacáo de 
dados através da automatizacáo dos processos de 
tratamento, uma grande quantidade de informacáo 
poderá ser gerada de forma sistemática e objetiva. 
Neste contexto o presente trabalho tem como objetivo 
a otimizacáo do processo de mapeamento 
automático de áreas urbanizadas de imagens do 
satélite Landsat, as quais seráo utilizadas em estudos 
relativos ao novo processo de urbanizacáo 
denominado urbanizacáo dispersa, que se 
caracteriza pela ocorréncia de pequenos núcleos 
afastados do tecido urbano contínuo. Com isto 
pretende-se diminuir a subjetividade associada a 
interpretacáo manual de dados de sensoriamento 
remoto, acelerar o processo de obtencáo de 
informacóes atualizadas e subsidiar políticas 
públicas e atuacáo de diferentes profissionais que 
tem como objeto de estudo o espaco urbano. O 
espaco urbano analisado neste estudo é o do 
município de Sáo José dos Campos, estado de Sáo 
Paulo, Brasil. 


Fenómeno Urbano 


Dada a importáncia e porte das modificacóes que 
estáo ocorrendo quanto a forma de crescimento das 
cidades brasileiras estáo sendo desenvolvidos 
estudos que apresentam um diagnóstico das áreas 
urbanizadas. Para Small (2005) embora as áreas 
urbanas ocupem uma porcáo relativamente pequena 
da superfície terrestre sua extensáo, distribuicáo e 
evolucáo tem grande impacto no meio ambiente e na 
dinámica sócio- económica mundial. 

Segundo Reis (2006) que desenvolve estudos das 
mudancas que estáo ocorrendo no processo de 
urbanizacáo do Brasil estas ocorrem : com a 
formacáo de áreas de urbanizacáo dispersa, que se 
estendem por vasto território, separadas no espaco 
físico mas mantendo estreitos vínculos entre si, como 
elemento de um único sistema; com a regionalizacáo 
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do cotidiano das pessoas que adotam novos modos 
de vida devido a mobilidade; com novas formas de 
gestáo do espaco pelo surgimento de formas 
condominiais diversificadas como condomínios 
residenciais, grandes shoppings centers, complexos 
universitários, entre outras. Os problemas dessa 
urbanizacáo sáo enumerados em estudo desen- 
volvido por Hasse e Lathrop (2003) para elaboracáo 
de medidas espaciais para melhor identificar e 
quantificar as características da dispersáo urbana. 

As fragilidades ecossistémicas segundo Carmo e 
Hogan (2006) exigem que no planejamento e na 
implementacáo de políticas públicas, a questáo da 
sustentabilidade seja priorizada desde o início. 

Ainda de acordo com estes autores os instrumentos 
de controle ambiental revelam-se frágeis frente ás 
press0es políticas. 

De acordo Maktav (2005), o crescimento náo 
controlado da populacáo e a urbanizacáo rápida 
especialmente em países em desenvolvimento tém 
grande impacto no meio ambiente podendo acarretar 
mudancas irreversíveis. Segundo Bonduki (2005) as 
cidades tém condicóes de planejar seu futuro, 
reduzirem as desigualdades, promover um adequado 
uso do solo para evitar desastres ambientais. Para 
que isto, o poder público deverá intensificar a 
fiscalizacáo da aplicacáo das leis , para que náo se 
transformem em meras formalidades que náo saíram 
do papel. 


Classificacáo orientada a objeto 


Conforme Lee and Lathrop (2006) o aumento da 
preocupacáo com a urbanizacáo dispersa de muitas 
zonas metropolitanas nos EEUU como também do 
globo tem provocado maior desenvolvimento de 
pesquisas interessadas no mapeamento e 
caracterizacáo do uso do solo urbano e suburbano, 
utilizando-se técnicas automáticas de interpretacáo 
dos dados. O problema do uso de imagens de 
resolucáo moderada para estes mapeamentos é a 
presenca de pixels mistos, compostos de superfícies 
impermeáveis, gramas, árvores resultando em uma 
mistura de respostas espectrais individualizadas. Os 
sistemas tradicionais de classificacáo de associar um 
pixel a um único uso do solo ou cobertura tornam-se 
problemáticos em áreas urbanas espacialmente 
heterogéneas com alta proporcáo pixels mistos. No 
trabalho desenvolvido por eles, que trata do 
desenvolvimento de sistema de classificacáo 
automática para caracterizacáo de tipos de 
coberturas urbanas, este problema é abordado. 

Segundo Maktav (2005) em adicáo aos 
melhoramentos tecnológicos na obtencáo de 
imagens de satélite as análises das imagens também 
tiveram que ser melhoradas. As abordagens 
orientadas a objetos estáo se mostrando melhores 
em muitas aplicacóes. A classificacáo orientada ao 
objeto consiste no agrupamento de pixels a partir da 
segmentacáo da imagem. Á segmentacáo é o passo 
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preliminar que divide a imagem em objetos 
homogéneos e contíguos. Á acurácia da 
segmentacáo afeta diretamente o desempenho da 
classificacáo. Segundo Geneletti and Gorte (2003) as 
técnicas de classificacáo orientadas a objetos 
baseadas em segmentacáo das imagens estáo 
despertando interesse como método de produzir 
mapas que sáo diretamente utilizáveis em Sistemas 
de Informacáo Geográfica (GIS). 

Nos sistemas de análise orientada a objeto os objetos 
se relacionam por meio da definicáo da rede 
hierárquica (heranca dos atributos que descrevem a 
classe) e da rede semántica (estrutura lógica de 
relacáo entre as classes). Para a classificacáo é 
utilizado o maior grau de pertinéncia dos objetos a 
determinada classe por meio da lógica fuzzy. 
Bruzzone and Carlin (2006) em trabalho desen- 
volvido para a criacáo de um sistema de classificacáo 
baseado em contexto multi-níveis para imagens de 
alta resolucáo espacial coloca que a classificacáo 
orientada a objeto e/ou estratégias multi-níveis/ multi- 
escala é baseada no sistema de comportamento 
humano na interpretacáo visual de dados. A 
racionalidade desta abordagem é que cada imagem é 
feita de objetos de diferentes tamanhos e formas 
interrelacionados. 

Assim cada objeto pode ser modelado pela forma e 
por medidas topológicas que podem ser usadas e 
integradas com características espectrais para 
melhorar a precisáo da classificacáo. Os objetos sáo 


extraídos das imagens com técnicas de segmentacáo. 


Em estudo desenvolvido por Lee and Warner (2006) 
na comparacáo entre classificacáo orientada a objeto 
com aquela pixel a pixel, a vantagem da primeira foi 
significante por combinar informacáo espacial e 
espectral já no processo de segmentacáo no qual os 
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objetos sáo criados usando as variáveis espectrais e 
espaciais. Neste procedimento foi diferenciada a 
classe área urbana construída. 

Tullis e Jensen (2003) desenvolveram trabalho cujo 
propósito foi demonstrar e testar um método 
automático para extrair unidades residenciais de 
dados multiespectrais de satélite, pois estes podem 
ser adquiridos para grandes áreas com alta resolucáo 
temporal. Foram utilizadas imagens IKONOS para 
deteccáo de unidades residenciais usando imagem 
NDVI, como um canal, e as variáveis: forma, tamanho, 
e contexto. Yan et al. (2006) testaram o desempenho 
da classificacáo baseada em pixel (Máxima 
Verossimilhanca) e da orientada a objeto para oO 
mapeamento de classes de uso e cobertura do solo 
visando a definicáo de áreas potenciais de queima de 
carváo em território chinés. O mapa temático 
resultante da análise orientada a objeto apresentou 
uma precisáo total de 83,25%, maior que aquela 
obtida pela classificacáo pixel a pixel, de 46,48%. As 
áreas construídas/edificadas foram somente 
discriminadas no mapa temático resultante da 
classificacáo orientada a objeto. Ouros estudos como 
os de Mitri and Gitas (2004) e Chubey et al. (2006), 
mostram a vantagem da classificacáo orientada a 
objeto. 


MATERIAL E METODOLOGIA 


Neste estudo foram utilizadas as imagens do satélite 
Landsat, órbita/ponto 219/76, bandas 3, 4 e 5 de 
12/03/1985 e 03/09/1999, referentes ao município de 
Sáo José dos Campos, localizado no Estado de Sáo 
Paulo (Figura 1). 





FIGURA 1 - Área de Estudo 
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Inicialmente, foi calculado o NDVI (bandas 4 e 3 do 
TM Landsat) visando uma melhor separacáo entre as 
classes de vegetacáo e as áreas urbanizadas. No 
mapeamento dos núcleos urbanos aplicou-se a 
classificacáo orientada a objeto utilizando o Software 
eCognition 5.0 (DEFINIENS, 2006) na classificacáo 
das imagens Landsat e da imagem NDV!I gerada. A 
análise orientada a objeto disponível neste software 
se baseia em 2 passos: Segmentacáo Multiresolucáo 
e Classificacáo. 

Em seguida foi feita a segmentacáo Multiresolucáo, 
baseada nos parámetros de escala, cor ou forma e 
textura, utilizando as bandas 3, 4 e 5 e a imagem 
NDVI. Optou-se pelo valor 10, com peso maior para o 
parámetro cor. Posteriormente, foi realizada uma 
segmentacáo denominada Diferenca Espectral, na 
qual os objetos obtidos na segmentacáo anterior 
foram fundidos em funcáo dos valores médios de 
intensidade que estes apresentavam em cada banda, 
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considerando uma distáncia espectral 10 (valor 
definido empiricamente). Em funcáo do resultado 
obtido repetiuse o procedimento de segmentacáo por 
Diferenca Espectral, com uma distáncia espectral de 
20, visando a melhor separabilidade do objeto de 
interesse, no caso, a mancha urbana. 

Esta abordagem de segmentacáo permite que se 
escolha o layer que melhor discrimina uma 
determinada classe. Os objetos obtidos na 
segmentacáo podem se relacionar por meio de uma 
relacáo semántica, hierárquica, ou ambas. Para se 
estabelecer estes relacionamentos entre as classes 
definidas para este estudo foram criadas duas redes, 
as redes semántica e hierárquica (Figura 2). As redes 
foram criadas em funcáo dos objetivos aqui propostos 
e das características espectrais e espaciais das 
imagens utilizadas. Os objetos de cada classe foram 
classificados de acordo com seus atributos de 
relacáo semántica, hierárquica, ou ambas. 
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FIGURA 2- Redes Semántica (a) e Hierárquica (b) 


A rede semántica mostra a relacáo entre as classes 
de interesse definidas no projeto, por exemplo: em 
áreas urbanizadas há vegetacáo rasteira e arbórea, 
corpos d 'água e solo exposto. 

Atributos relacionais podem ser criados para a 
caracterizacáo de uma determinada classe. A rede 
hierárquica define as herancas de atributos “físicos” 
entre as classes, tais como características espectrais, 
espaciais de forma, textura, etc. 

A próxima etapa é a determinacáo dos atributos que 
descrevem as classes de interesse. Foram utilizados 
aqueles atributos que apresentaram melhores 
resultados na discriminacáo das classes. 

Cada classe foi determinada por atributos de forma, 
cor, textura e contexto que, isoladamente ou em 
conjunto, melhor as caracterizam, permitindo um 
melhor resultado na classificacáo final. No software 
utilizado, eCognition 5.0, existem diversos atributos 
pré-definidos, podendo-se criar um número indefinido 
de atributos relacionais ou aritméticos a partir das 
características do objeto que o intérprete julgar 
relevantes. 


Todos os atributos, disponíveis no eCognition, ou 
criados a partir da necessidade de classificacáo e 
características dos dados, foram espacializados para 
que fossem selecionados aqueles que melhor 
discriminam os objetos de uma determinada classe 
emrelacáo a toda a área de estudo. 

Com base nessa espacializacáo foram selecionados 
os atributos que melhor caracterizavam as classes 
por meio da análise visual. Na imagem Landsat 
segmentada foram selecionados os polígonos 
amostrais referentes a cada classe analisada. 

A partir das amostras foram gerados gráficos 
utilizando os valores dos atributos de cada classe. 
Comisso foram obtidas as funcóes de pertinéncia dos 
atributos descritores de cada classe. Esta análise 
auxilia tanto no processo de selecáo dos próprios 
atributos e na definicáo do range de cada um deles, 
quanto no refinamento da selecáo das amostras de 
cada classe. Nesta etapa, os valores obtidos para as 
amostras das diferentes classes foram comparados 
visualmente nos gráficos gerados e por meio de 
valores de overlap e desvio padráo. Para a classe 
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área urbanizada, de interesse deste trabalho, foram 
utilizados os atributos de nível de cinza na banda 3 e 
na imagem NDVI e homogeneidade de textura. 

O atributo de textura foi utilizado para a 
separabilidade entre as classes área urbanizada e 
solo exposto, devido a confusáo que ocorre entre 
estas duas classes. Dentre os atributos de textura 
disponíveis foi selecionada a feicáo homogeneidade, 
já que as áreas urbanizadas apresentam-se com um 
padráo rugoso em relacáo as áreas com solo exposto. 
Como etapa final do processo de mapeamento das 
áreas urbanizadas foi feita a classificacáo a partir dos 
atributos e seus intervalos de valores, definidos por 
meio das amostras selecionadas. O mapeamento por 
meio do sistema de classificacáo orientada a objeto 
apresentado foi realizado para as imagens de 1985 e 
1999. 

A avaliacáo da precisáo da classificacáo, para a área 
de estudo foi feita, para o ano de 1985, a partir da 
comparacáo do mapa obtido pela interpretacáo 
automática com o mapeamento do município obtido 
por meio da interpretacáo manual de fotografias 
aéreas na escala 1:25.000 (Sáo José dos Campos, 
Prefeitura Municipal, 1985), em funcáo da 
disponibilidade de um mapeamento obtido por meio 
da interpretacáo manual de fotografias aéreas na 
escala 1:25.000. Para este cálculo as classes 
analisadas foram agrupadas em áreas urbanizadas e 
náo urbanizadas. 

Os resultados da classificacáo do ano de 1999 foram 
avaliados tendo como base de referéncia o mosaico 
aerofotogramétrico digital do ano 2000 (Sáo José dos 
Campos, Prefeitura Municipal) da área urbanizada de 
Sáo José dos Campos. Esta área corresponde a 
maior parte do perímetro urbano e de expansáo 
urbana que ocupa a regiáo sul do Município. 

Para a avaliacáo dos dados foram utilizadas 
quadrículas amostrais de 1Km por 1Km superpostas 
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ao mosaico aerofotogramétrico. Os erros de omissáo 
e inclusáo foram identificados a partir da comparacáo 
visual entre os resultados da classificacáo digital e o 
mosaico. O critério utilizado na avaliacáo das 
quadrículas foi o de considerar o erro dominante em 
situacóes em que foram observados os dois tipos de 
erros. 


RESULTADOS 


Durante o processo de classificacáo verificou-se que 
na classe Áreas Urbanizadas foram incluídas áreas 
de Solo Exposto. Isto ocorre devido ás semelhancas 
espectrais destas duas classes, pois mesmo com a 
utilizacáo do atributo de textura, que reduz, mas náo 
elimina totalmente a confusáo entre elas, nas áreas 
em que o relevo é acentuado, o sombreamento nas 
encostas fornece valores de textura semelhantes 
para ambas. Neste caso, deve-se ter cautela, pois 
áreas de transicáo de uso com solo exposto, como 
loteamentos novos, áreas industriais e shopping 
centers em implantacáo correspondem a classe de 
interesse, ou seja, urbana. 

Estes erros poderiam ser reduzidos utilizando-se 
dados auxiliares no processo de classificacáo como, 
por exemplo, uma malha viária completa. Neste caso, 
a proximidade em relacáo as estradas poderia ser 
utilizada como um atributo de restricáo espacial na 
classificacáo A classe área urbanizada foi separada 
da classe vegetacáo por meio do atributo NDVI, 
inclusive no interior das áreas urbanizadas (mata 
ciliar, parques, bosques, pragas). 

Os resultados finais para a classe Area Urbanizada, 
obtidos nas classificacóes sobre as imagens dos 
anos de 1985 e 1999, podem ser verificados na 
Figura 3. 
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FIGURA 3 Classificacoes obtidas por análise orientada a objeto, realizada no eCognition, sobre imagens 
Landsat 5/TM de 1985 (a) e 1999 (b) 
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Pela análise dos resultados apresentados na Figura 3 
verificou-se que as áreas agrícolas com solo 
preparado para plantio náo foram classificadas como 
urbanas como ocorre nas classificacóes pixel a pixel 
em que é considerada apenas a característica 
espectral dos alvos. Outro aspecto observado é que 
muitas áreas de pastagem em relevo colinoso com 
elevada porcentagem de solo exposto gerando 
resposta espectral semelhante ás áreas urbanizadas, 
embora de textura diferente, náo foram incluídas 
como classe urbana. 

Os loteamentos com baixa densidade de construcáo, 
como aqueles que apresentavam um padráo de 
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loteamento de chácaras, náo puderam ser 
classificados como áreas urbanizadas devido aos 
altos índices de vegetacáo na imagem de 1985 
(Figura 3a). Isto pode ser resolvido dentro da própria 
rede hierárquica, subdividindo-se a área urbanizada 
em consolidada e náo-consolidada. 

A Tabela 1 mostra o resultado da comparacáo entre a 
classificacáo realizada no eCognition para o ano de 
1985 e o Mapeamento de Uso da Terra de Sáo José 
dos Campos quanto á classe Área Urbanizada. As 
demais classes obtidas foram agrupadas em Náo- 
Urbanizada. 


TABELA 1 - Matriz de confusáo da Classificacáo Orientada a Objeto da imagem Landsat de 1985 


Classes de Uso da 
yy Un rbanizada ES 





A classificacáo é considerada como muito boa, de 
acordo com a tabela de conceitos de eficiéncia dos 
índices Kappa. Os erros de omissáo e inclusáo 
representam, além dos erros de classificacáo 
propriamente ditos, a diferenca na identificacáo dos 
usos nos sistemas manuais e automáticos de 
classificacáo. Essa dificuldade se dá devido a 
diversidade dos alvos urbanos que ocorre em um 
único uso funcional do solo urbano. No caso da 
interpretacáo de dados de alta resolucáo osdiferentes 
usos sáo mapeados considerando-se todos os 
elementos que podem constituí-lo como áreas 
construídas, vegetacáo e água. 

Quanto ao Mapa de Uso da Terra do Município de Sáo 
José dos Campos (Sáo José dos Campos, Prefeitura, 
1985), utilizado como referéncia, os limites das 
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diferentes classes de uso do solo foram determinados 
a partir da interpretacáo de fotografias aéreas, 
representando, na maioria dos casos, os limites dos 
lotes. No método automático de classificacáo aqui 
avaliado, estas propriedades sáo representadas 
pelos limites de suas coberturas naturais do solo e 
áreas impermeabilizadas, já que o uso do solo náo 
pode ser inferido, principalmente devido a resolucáo 
espacial das imagens de satélite aqui utilizadas. 

A Figura 4 mostra a distribuicáo espacial dos erros de 
omissáo e inclusáo resultantes da classificacáo da 
imagem do ano 1999. De um total de 423 quadrículas 
analisadas, 329 tiveram concordáncia com o mosaico 
de referéncia utilizado; 26 apresentaram erros de 
omissáo, enquanto 68 apresentaram erros de 
inclusáo. 
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FIGURA 4 - Distribuicao espacial dos erros de omissáo e inclusáo obtidos na classificacáo da imagem de 1999 
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Os erros de omissáo ocorreram em menor número 
em relacáo aos erros de inclusáo e encontramse 
associados: a áreas edificadas junto a estradas 
vicinais com arruamento formado apenas por uma ou 
duas vias, constituindo-se em aglomeracóes de 
desenvolvimento linear junto a estradas que 
funcionam como ponto inicial e a partir das quais se 
desenvolvem transversalmente; a áreas em 
consolidacáo com baixa densidade de residéncias. 
Os erros de inclusáo que provocaram super- 
estimativa das áreas urbanizadas foram ocasionados 
pela classificacáo de áreas de pastagem, na maior 
parte dos casos, ou campo antrópico com edificacdes 
esparsas. 

Grandes áreas em extensáo territorial adjacentes ás 
edificacóes isoladas, principalmente no norte do 
município foram classificadas como urbana. Elas 
correspondem a áreas ocupadas por vegetacáo de 
gramíneas podendo ser áreas de pastagem ou 
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campos antrópicos e de relevo colinoso com 
vertentes iluminadas. 

Embora ocorra a domináncia de áreas náo 
urbanizadas estas áreas apresentam ocupacáo 
rarefeita em edificacóes isoladas. Junto a área 
urbanizada consolidada principal verificaram-se 
problemas de classificacáo relativos a áreas em que 
há presenca de solo exposto devido a erosáo. 

Na franja rural urbana houve erro relacionado a 
classificacáo de áreas de granjas. O material utilizado 
na cobertura dos galp0es usados para a confinacáo 
das aves apresentou a mesma resposta espectral e 
forma regular que as áreas urbanizadas. 

Na regiáo sul do município a rotatória do 
entroncamento de rodovias intermunicipais foi 
classificada como urbana caracterizando assim erro 
de inclusáo. A Tabela 2 mostra os erros de omissáo e 
inclusáo relativos á classificacáo do ano 1999. 


TABELA 2 - Matriz de confusáo da Classificacáo Orientada a Objeto da i¡magem Landsat de 1999 








A classificacáo é considerada boa, de acordo com a 
tabela de conceitos de eficiéncia dos índices Kappa 
0.55. O índice inferior deste em relacáo áquele 
referente ao obtido na classificacáo da imagem de 
1985, pode estar relacionado ao procedimento 
utilizado, bem como aos dados de referéncia. Um 
outro fator que pode explicar esse resultado, é que 
para 1985 foi considerado todo o município e, 
consequentemente, há uma ocorréncia maior de 
áreas náo urbanizadas classificadas corretamente. 


CONCLUSOES 


Os resultados mostraram um bom desempenho das 
classificacdes realizadas ao separar os objetos que 
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